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It is widely accepted that moduli in the mass range lOeV - lO^GeV which start to oscillate with an 
amplitude of the order of the Planck scale either jeopardize successful predictions of nucleosynthesis 
or overdose the Universe. It is shown that the moduli problem can be relaxed by making use of 
parametric resonance. A new non-perturbative decay channel for moduli oscillations is discussed. 
This channel becomes effective when the oscillating field results in a net negative mass term for the 
decay products. This scenario allows for the decay of the moduli much before nucleosynthesis and, 
therefore, leads to a complete solution of the cosmological moduli problem. 
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I. INTRODUCTION 

In many theories beyond the Standard Model of 
particle physics, in particular in supergravity and 
string theories, there are many scalar and fermionic 
fields with masses smaller or equal to the elec- 
troweak scale and gravitational strength couplings 
to ordinary matter. Such fields are called mod- 
uli fields and behave as non-relativistic matter at 
late time. Since they decay late because of their 
weak interactions, they lead to the so-called cos- 
mological moduli problem [1, 2, 3, 4]. Below, inde- 
pendently of their origin, fields having only Planck 
scale couplings and weak scale mass are denoted 
collectively as moduli fields. 

In the case of moduli in supersymmetric theo- 
ries, the mass of these fields (j) is generated 
during supersymmetry breaking. A wide variety 
of these scenarios predict masses of the moduli in 
the dangerous range, ~ lOeV — W^GeV. Ac- 
cording to our ideas of early Universe cosmology, 
we expect moduli to be produced in great abun- 
dance in the early universe. In the context of Big 
Bang cosmology, both scalar and fermionic moduli 
particles will be part of the initial thermal bath of 
particles of the very early universe. Even assuming 
that the moduli particles are not part of the initial 
thermal bath (for example in the context of infia- 
tionary cosmology) it is hard to avoid the presence 
of excited moduli fields at late times. For exam- 
ple, in the case of scalar moduli, since the moduli 
are massless before supersymmetry breaking, there 
is no reason that the moduli field values before 
supersymmetry breaking coincide with the values 



which turn into the minima of the potential after 
supersymmetry breaking [5]. The offset will lead 
to moduli fields which oscillate about their poten- 
tial minima. An offset of a scalar modulus can also 
be produced by quantum fiuctuations in the early 
phases of infiation, as follows from the computa- 
tion of the coincident point two point function of 
a low-mass scalar field during inflation [6, 7]. The 
excessive production of moduli is predicted during 
the preheating stage of inflationary cosmology in 
a wide variety of models [8]. A further source of 
moduli particles is gravitational particle produc- 
tion between the end of inflation and the time of 
nucleosynthesis [9, 10]. 

Due to their weak interactions, the decay of the 
moduli fields is slow. Widely used estimates based 
on dimensional analysis for the perturbative decay 
rate T give 



The presence of excited moduli fields at late times 
is dangerous since the presence of the extra mod- 
uli field energy during the time of nucleosynthesis 
could spoil the success of the standard Big Bang 
nucleosynthesis scenario [11]. This danger is acute 
in particular for the heavier moduli fields. Both 
lighter moduli fields and heavier fields which do 
not decay before the time of equal matter and ra- 
diation threaten to overdose the Universe at that 
time (and thus also today if they do not decay be- 
tween the time of equal matter and radiation and 
the present time). 

To be more specific, it is the fact that 
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the interactions of moduli fields come from 
non-rcnormalizablc terms in the supersymmetry 
(SUSY) potential that leads to the moduli prob- 
lem. For example, the following term in the Kahler 
potential 

^K=^cl^\h<l>'^<l>, (2) 
p 

results in a contribution to the square mass for (j) of 
the form p/M^ [12]. By ((> we denote a canonically 
normalized modulus field (with bare mass m,^), by 
Mp - the Planck mass, by 0/ - a field which dom- 
inates the energy density of the universe, the in- 
flaton, and by p - the energy density contained in 
the infiaton. Since p = SH'^Mp, the effective mass 
nieff becomes 

ml^.j = cH^ + ml, (3) 

where c is a constant. If c 1, then inflation 
drives the moduli to the minima of their high tem- 
perature effective potential. However, typically the 
high temperature minima are offset from the zero 
temperature minima of the moduli potentials by 
a value which has Planck order of magnitude. If 
c <C 1 then H > m^f f , and then quantum fluctu- 
ations during inflation will also excite the fleld (j) 
to a value of the order of the Planck mass. Dur- 
ing reheating the energy density of the inflaton 
stops dominating the Universe and the effective 
mass of (j) relaxes to m^. After reheating, the 
Hubble constant decreases in the radiation dom- 
inated phase. Once H ^ m.^ the condition for 
slow rolling of </> (V^"(^) < H^) is no longer satis- 
fied, and at that point the field (p starts to oscil- 
late around its low temperature minimum which 
we take to be zero. The energy density of (j) de- 
creases like that of non-relativistic particles (i.e. 
proportional to a{t)~^, where a{t) is the cosmolog- 
ical scale factor) , whereas the energy density of the 
radiation dominated universe falls as a{t)^*. Thus, 
the modulus field (p may come to dominate the en- 
ergy density of the universe, or at least contribute 
too much during the period of nucleosynthesis, un- 
less the field decays early. The problem is that the 
modulus fields are coupled only gravitationally to 
themselves and other particles and thus have very 
small perturbative decay rates. 

There has been some previous work to try to 
mitigate the cosmological moduli problem. En- 



tropy production at late times will dilute the mod- 
uli density. For example, weak scale inflation [13] 
or thermal inflation [14] could sufficiently dilute 
unwanted moduli (see e.g. [15] for a detailed study 
of the potential of thermal inflation to solve the 
moduli problem for masses in the above-mentioned 
dangerous range predicted in models with hidden 
sector and gauge-mediated supersymmetry break- 
ing). For certain ranges of parameters of a gauge- 
mediated supersymmetry breaking model, the os- 
cillations of the modulus field itself might suffi- 
ciently dilute the string moduli density [16]. The 
decay of an unstable domain wall network [17] is 
another way to generate entropy and dilute the 
moduli density at late times. A common danger of 
these approaches is that the baryon density might 
also be diluted to an unacccptably low value. An- 
other approach is to invoke effects which give the 
moduli fields a contribution to the square mass of 
the order of which would allow them to roll 
down their potential during inflation [18, 19] and 
prevent them from acquiring a large expectation 
value during inflation (see also [20]). However, this 
solution does not work if the low temperature min- 
imum of the moduli potential does not coincide [5] 
with its high-temperature ground state (new sym- 
metries which could force the two states to be the 
same were analyzed in [12]). A recent proposal to 
solve the moduli problem is moduli trapping at en- 
hanced symmetry points [21, 22]. In terms of the 
use of parametric resonance instabilities, our work 
has similarities with that of [21, 22]. However, in 
contrast to these works, in our study the focus is 
on the traditional moduli problem as formulated 
in [1, 2, 3, 4]. 

In this paper, we propose a way to solve the cos- 
mological problems of scalar moduli fields which 
requires no external mechanism for the dilution of 
moduli. Instead, it makes use of non-perturbative 
decay channels. Non-perturbative decays have 
been shown to completely change the scenario of 
reheating in inflationary cosmology. In particular, 
the decay of the inflaton field by a parametric reso- 
nance instability has been shown to be very impor- 
tant [23, 24, 25, 26]. In certain models, a tachyonic 
instability renders the inflaton decay even more ef- 
ficient [27]. 

We consider a toy model which potentially 
gives rise to a cosmological moduli problem. In 
the framework of this model we investigate non- 
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perturbative decay channels (decays into particles 
whose lifetime is much shorter than that of the 
moduli fields). We study the decay of the oscil- 
lating modulus field via parametric resonance, and 
propose a new tachyonic decay channel. We quan- 
tify the conditions on the parameters of the model 
for which the decay channels are eff'ective. Note 
that the decay channels work for an initial field 
amplitude up to the order of the Planck scale. No 
external mechanism for diluting the moduli density 
is required. 



II. THE MODEL 



In the above, the mass mi sets the scale of the 
interaction between (f) and x- It is reasonable to 
assume that this scale is much less than Mp. 

In analogy to the situation encountered in the 
study of the decay of the inflaton, if we want to 
study the decay of the modulus field, it is crucial to 
focus on the equation of motion for matter fields x 
which the modulus field couples to. In the presence 
of the oscillating modulus field, this equation (in 
an expanding space-time) is 

Xk + SHick + -j^Ht) sin(m^t) j Xk 

= 0, (8) 



If the moduli problem arises as a consequence 
of supersymmetry breaking, the moduli potential 
takes the form (see e.g. [13]) 



V($) = ml/.M^gmMp) 



(4) 



where m3/2 is the gravitino mass and Q is some 
function. We will assume that the modulus field 
couples to some matter field x which we treat as 
a scalar field (following the analyses in the study 
of fermionic preheating [28] it could also be taken 
to be a fermionic field). Making use of (4), and of 
the fact that the dimension five and six operators 
are suppressed by Mp and by Mp, respectively, we 
can construct a typical potential for the modulus 
(j) and the matter scalar field x 



V = 



I 2/2,1 22 

+Xi(l)\^ + \2X^ + A3<; 



I'Tl? 2 

2Ai^^^ 



(5) 



where the coupling constants Ai, A2, A3 are small 
enough such that for \(f)\, \x\ < Mp the low energy 
efii'ective potential for cj) takes form 



Vi = -m; 



, 1 2 2 , l™r / 2 
+ 2™>^^ +2M/^ 



(6) 



Thus, the toy model we consider is described by 
the following Lagrangian: 



d''<f>d^<j) + d^'xd^x 
1 



(7) 



--m 
2 



1 2 2 ; 2 

2"^-^ -2M-p^'' ■ 



where $(t) is the amplitude of (j). The ampli- 
tude $ decreases as a consequence of the expan- 
sion of space. In the above, k denotes the comov- 
ing wavenumber and a{t) is the cosmological scale 
factor. 

In a first step, we will put the above equation 
into the form of the well-known Mathieu equa- 
tion. To absorb the expansion of space, we define 
a rescaled field via 



Vk = a^^'^Xk' 



(9) 



Then, the equation of motion for rik becomes: 



rik+{-^ + ml + jf^t) sin(m^t) - A ) 77^ = , 



where 



(10) 



^=l"' + l-a = ^(l/2 + (p-l)/p)i^^ (11) 

where for the second equality we have assumed 
that a{t) cx fP. 

It is convenient to introduce a dimensionless 
time variable via 



1 TT 

z = -m^t + - 



(12) 



The diflferentiation with respect to z will be de- 
noted by a prime. In this case, the above equation 
(10) takes the form 



r]',^ + {Ak-2qcos2z)r]k = 0, 



(13) 
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where 



I u 4 ± 

mlo^ ml 



4m^^A 



(14) 
(15) 



If it is justified to neglect the expansion of space 
compared to the rate of the processes which will be 
discussed in the following, then A vanishes, and 

the equation (13) takes the form of the Mathicu 
equation. Note that fc/a is the physical wavenum- 
ber. 



III. PARAMETRIC RESONANCE 
INSTABILITY 

The modulus field (p is frozen until the Hubble 
parameter H drops to a value comparable to m^. 
Then, begins to oscillate about = with a fre- 
quency m0, its amplitude <I>(t) being damped by 
the expansion of the spatial background and by the 
energy loss of cp to other fields. The second effect 
is a back-reaction effect which we will neglect. The 
condition required that the oscillation of (f> begins 
before the time of nucleosynthesis is that the Hub- 
ble damping term in the equation of motion for (p 
becomes smaller than the force term V (</>) driving 
the oscillations. It yields 



HiTNs) < m^, 



(16) 



(where Tjvs is the temperature at which nucleosyn- 
thesis takes place) a condition which is satisfied for 
all masses in the dangerous range. Once the mod- 
ulus field starts to oscillate, resonant excitation of 
all fields coupled to (p is possible, in particular the 
excitation of rj. 

The first instability we will study is the paramet- 
ric resonance instability, first applied to the decay 
of the inflaton field in [23]. There are two types 
of resonance [24], broad parametric resonance and 
narrow parametric resonance. The condition for 
broad resonance is g > 1. where q is the parameter 
appearing in equation (13). This condition is sat- 
isfied for large values of the amplitude namely 
for 



(17) 



Evaluating this condition at the time when pertur- 

bative decay of the modulus field sets in, i.e. when 
r = if, and using for formula for F applicable to 
our toy model (given below in (29)), we find that 
unless 



<C m/( 



mi si/2 



10M„ 



(18) 



broad parametric resonance can relax but not solve 
the moduli problem without an additional decay 
channel being present. Both narrow resonance 
(disciisscd below) and the tachyonic decay dis- 
cussed in the following section can provide the ad- 
ditional channel. 

For smaller values of $, we are in the domain of 
narrow parametric resonance. In this phase, the 
growth of rik is known [29, 30]: 



•nk 



(19) 



Only modes in narrow resonance bands are ampli- 
fied, and the first such band is centered at a value 
of k = km given by [29, 30] 



1 , kn 



^(1 



ml 



(20) 



from which it follows that the band does not exist 
unless < m,0/2. Other resonance bands occur 
for larger values of but are of higher order in 
perturbation theory and hence have a negligible 
effect. 

The first condition for resonance to be effective 
is that the typical time scale for the growth of r]k 
is shorter than the Hubble time, i.e. 



qm^ > H. 



(21) 



In addition, one should take into account the 
change of the momentum as a result of the back- 
ground expansion [26]. In a time interval St, as- 
suming TO^ < m^/A, the change in the physical 
momentum p = k/aif) corresponding to the mid- 
dle of the lowest resonance band is 



Ap = pHM 



HM. 



The width of the resonance band is 
Ap=^. 



(22) 



(23) 
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Thus, p remains in the resonance band during a 
time interval 



At ~ qH- 



(24) 



To justify neglecting the expansion of space, we 
must require that the exponent in the growth fac- 
tor (19) is at least 1 during this time interval. This 
leads to a more severe constraint on q: 



> H. 



(25) 



Inserting the value of q from (15), we find that 
narrow resonance is efficient provided 



$(t) > $e(t) = 



HMpm 



3/2 



rrij 



(26) 



Since decreases as t ^1"^ whereas in a 

radiation-dominated phase $(t) decreases only as 
narrow resonance decay channel eventu- 
ally shuts off. In order that the moduli field does 
not dominate the energy density at the time of the 
shutoff at temperature T, the following condition 
needs to be satisfied: 



(27) 



which, inserting the expression (27) for $c, turns 
into the condition 



(28) 



For narrow parametric resonance to solve the 
cosmological moduli problem, one needs to check 
that at the time of moduli decay (which occurs 
when r ~ iJ), the condition (28) still holds. The 
decay rate of ^ for our Lagrangian (7) is given by 



327rM2m^ " 



(29) 



GeV, the problem is solved for values of satis- 
fying 

IQ-'^GeV <m^< IQ-'^GeV . (31) 

In the context of our model (for mj = lO^GeV), 
moduli with masses smaller than 10~^GeV decay 
before time of nucleosynthesis and thus do not 
cause the problem. Note that this scaling of the 
decay rate with the mass of the decaying particle 
is going against the intuition that lighter moduli 
shoiild decay later than heavier one. This curious 
aspect of our toy model decreases the potentially 
dangerous mass range, and this realization might 
be useful in some concrete models suffering from a 
moduli problem. 

It appears at this point of our study that the 
period of narrow parametric resonance has the po- 
tential of solving the moduli problem for values 
of and to/ which satisfy the relation given by 
(30). One issue which we have not taken into ac- 
count is the fact that late moduli decay may pro- 
vide a large source of non-thermal photons which 
could distort the black-body nature of the CMB. 
The constraints resulting from this effect must be 
studied in any concrete model with late-decaying 
moduli fields. 

The previous analysis has missed a second im- 
portant condition for the efficiency of narrow para- 
metric resonance. It is not sufficient that the 
modes rjk increase with a rate faster than H. Since 
the resonance occurs only in narrow bands [29, 30], 
it is important to check that the rate of energy in- 
crease integrated over all modes of rj be larger than 
the decrease in the energy density of (f) taking into 
account the expansion of space alone. Otherwise, 
the energy stored in the moduli field would still 
scale as matter in spite of the exponential increase 
in the occupation number of certain field modes. 
This condition reads 



(32) 



Thus, the solution of the cosmological moduli 
problem requires: 



mi si/3 



(30) 



Inserting the temperature of nucleosynthesis to get 
lower bound from (29) on the potentially denger- 
ous mass range, and using the value mi = 10^ 



The rate of increase in the energy density of 
rj can be estimated by considering the increase in 
the amplitude of all modes of rj in the lowest in- 
stability band. This band is located at /c ^ m^ 
and its width is given by qm^. Since the rate of 
increase (from (19)) is qm^ and since the initial 
mode (vacuum) energy is about k, we obtain 



m%q^- 



(33) 
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Hence, the condition (32) for efficiency of the res- 
onance process becomes 



< mi 



mj mi 



(34) 



which is to be evaluated at the temperatures when 
the presence of moduU fields are dangerous for cos- 
mology. Using, as before, the value m/ — lO^GeV, 
and evaluating the above condition at the temper- 
ature of nucleosynthesis, we find that the condi- 
tion is satisfied as long as the mass m^ is smaller 
than about lO^'^GeV. The condition (34) becomes 
increasingly well satisfied at lower temperatures, 
and is no longer a concern at the time of recom- 
bination. Note that the condition (32) for the ef- 
ficiency of narrow resonance is a very conservative 
one. As long as the first condition (25) is satisfied, 
the mode amplitude will grow, and hence the vac- 
uum mode amplitude used in the above estimate 
should be replaced by the excited amplitude, thus 
relaxing the constraint by an exponential factor. 

The problem, however, is that the condition (34) 
conflicts up to coefficients of order unity with the 
condition that the perturbative decay rate is neg- 
ligible. Thus, in our toy model, and using the very 
conservative form of our conditions for the effec- 
tiveness of the resonance process, the narrow reso- 
nance decay channel is only effective near the time 
when the perturbative decay is also become im- 
portant. This result, however, is a consequence of 
the particular scaling of F with m^, T oc m^^. fn 
models where moduli decay to fermions, F a m^, 
and, therefore, we expect those models do not suf- 
fer from this specific problem. 

The above discussions neglected the expansion of 
the universe. Taking into account this expansion 
changes the Mathieu equation into a more general 
equation of Floquet type, and leads to a stochastic 
nature of the resonance process [26]. However, the 
property that the number of particles is growing 
exponentially at a rate given by (19) is preserved. 



IV. TACHYONIC DECAY OF THE 
OSCILLATING MODULUS FIELD 

In the case of the decay of the inflaton field at 
the end of the period of inflation, it is known [27] 
that for certain models there is a tachyonic insta- 
bility channel which is more efficient than para- 



metric resonance. In this section, we will study a 

similar process for moduli decay. 

Let us return to the basic equation (13), with 
the values of the parameters Ak and q given by 
(14) and (15), respectively. We immediately see 
that for large values of $, the effective term 
in the equation will be negative for part of the 
oscillation period of (j). This tachyonic instability 
occurs provided 



> 



m: 



m 



(35) 



The minimal value for which the tachyonic decay 
channel is open is given by setting the two sides in 
(35) equal and will be denoted by 

The condition under which the tachyonic decay 
channel can solve the moduli problem is then given 
by the condition 



^2 ^2 



$1 < 



(36) 



where T is the temperature corresponding to the 
period one is interesting in. Evaluating (36) at the 
time of nucleosynthesis, we obtain 



Tns mj 

Mr, ml 



(37) 



The upper mass bound on m^ for which the above 
tachyonic decay is effective thus depends sensi- 
tively on the ratio of m/ and . Unless the latter 
mass is much smaller than the former, the tachy- 
onic decay channel cannot reduce the amplitude of 
moduli oscillations to a level consistent with the 
observational constraints. 

The range of values of m^ for which the two 
decay channels - narrow parametric resonance (ne- 
glecting for a moment the issue that in our model it 
starts to be efficient together with the perturbative 
decay) and tachyonic decay - are open depends on 
the values of the masses. While if m^ > to^, the 
only allowed channel is the tachyonic one, the nar- 
row parametric resonance works for a wider range 
of masses when ~ m^. The later can be 
seen by inserting m^ = m^ into (37) and compar- 
ing with (28). However, as the value of m^ is re- 
duced, the range of masses for which the tachyonic 
decay channel is open grows and begins to dom- 
inate over that of narrow parametric resonance. 
Tachyonic decay can be made to work even for 
mi ~ m^. 
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For values of for which both decay channels 
arc open, the tachyonic decay is more efficient for 
two reasons. First, it leads to the excitation of all 
long wavelength modes, and not just to modes in 
a narrow resonance band. Modes with 

fcp < = kcrit (38) 

are excited. For such modes, the value of rjk in- 
creases with a maximal rate given by 

% ~ exp{y/qz), (39) 

which for (7 < 1 is a larger rate than that which oc- 
curs for modes in the resonance band during nar- 
row resonance (see (19)). This is the second rea- 
son for the larger efficiency of the tachyonic decay 
channel. 

Let us estimate the energy density stored in 
the quanta produced during the tachyonic decay 
process. The phase space of modes which are ex- 
cited tachyonically is of the order k^^.^^ . Each mode 
grows with a rate which varies in time, the maxi- 
mal rate being given by (39), and the growth oc- 
curs for approximately half the oscillation period 
(the period during which the effective square mass 
is negatiye). The mean growth rate is given by 
l/\/2 of the maximal rate. Thus 

Priit) ~ mij^exp{—mi^/¥/Mpt). (40) 

The prefactor in front of the exponential factor is 

much larger than the corresponding factor in the 
case of the narrow resonance decay process. 

V. DISCUSSION AND CONCLUSIONS 

In this paper we have studied two non- 
perturbative decay processes which can substan- 
tially dilute the density of dangerous moduli fields. 
Both processes occur during the phase in which the 
moduli fields oscillate about their ground states. 
The first non-pcrturbativc process is parametric 
resonant excitation of fields coupled to the modu- 
lus field. In this paper we modeled such fields as 
a scalar field x coupled to the modulus field (j> via 
dimension five operators suppressed by the Planck 
mass. The second decay process is tachyonic de- 
cay which makes use of the fact that the effective 



square mass of x in the presence of the oscillating 
(p field is negative for part of the oscillation pe- 
riod. These decay processes are analogous to the 
parametric resonant decay of the inflaton during 
reheating [23], and to tachyonic preheating [27], 
respectively. 

We have established the conditions under which 
either of the two decay processes can solve the 
moduli problem, i.e. reduce the energy density 
of the modulus field to values consistent with big 
bang nucleosynthesis. It appears that narrow para- 
metric resonance has the potential to solve the 
modulus problem but, for the toy model considered 
here, and using very conservative conditions for the 
effectiveness of the resonance process, it happens 
only at the time when the pertubativc decay rate 
also becomes important. Undoubtedly, tachyonic 
decay successfully solves the problem given masses 
of decay products which are much smaller than 
the scale of interaction - to/. Moreover, the tachy- 
onic decay channel allows for excitation of particles 
with masses heavier than that of the decaying par- 
ticle. Depending on the values of the other masses 
in the Lagrangian, either of the two decay pro- 
cesses can be open for a wider range of masses m^. 
For values of for which both decay channels are 
open, the tachyonic decay is much more efficient, 
as is true in the case of the decay of the inflaton. 
Note that the presence of an interaction term in 
the Lagrangian linear in (j) was important in order 
to obtain the tachyonic decay. 

It will be of great interest to study the applica- 
bility of these decay channels to concrete models 
with moduli flelds. This work is left for future re- 
search. 

The form of the potential in our toy model La- 
grangian inevitably suggests inflation at low scales. 
This natural source of inflation does not only di- 
lute heavier relics but could also mitigate the flat- 
ness, horizon and entropy problems. It needs to be 
studied whether this type of models can provide a 
successful reheating mechanism, if the infiaton has 
only gravitationally suppressed interactions. Once 
again, non-perturbative instabilities like those used 
in preheating [23, 24, 25, 26, 27] are likely to be 
successful. If the modulus field comes to domi- 
nate the energy density of the universe for some 
period (without necessarily leading to infiation), it 
can provide a candidate for the curvaton (see e.g. 
[31] for an extensive discussion of moduli fields as 
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candidates for the curvaton). 
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